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Abstract: 

Coupled asymmetric double well (acfr 2 — bcj) 5 + ctfi 4 ) one-dimensional potentials arise in the context of first 
order phase transitions both in condensed matter physics and field theory. Here we provide an exhaustive 
set of exact periodic solutions of such a coupled asymmetric model in terms of elliptic functions (domain 
wall arrays) and obtain single domain wall solutions in specific limits. We also calculate the energy and 
interaction between solitons for various solutions. Both topological (kink-like at T = T c ) and nontopological 
(pulse-like for T ^ T c ) domain wall solutions are obtained. We relate some of these solutions to domain 
walls in hydrogen bonded materials and also in the field theory context. As a byproduct, we also obtain a 
new one parameter family of kink solutions of the uncoupled asymmetric double well model. 
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1 Introduction 

First order phase transitions are typically modeled by a triple well (<f> 6 ) free energy. However, if a third 
order term becomes symmetry allowed, one need not go to the sixth order term for the transition to 
be of first order. Instead an asymmetric double well (or (ft 2 -eft 3 -(ft 4 ) free energy is sufficient to drive the 
transition (with or without an external field ^Hj). This situation occurs in body-centered cubic (bec) 
to face-centered cubic (fee) reconstructive structural phase transitions in crystals, the w-phase transition 
in various elements and alloys |2], isotropic to nematic phase transition in liquid crystals [21 and in the 
hydrogen chains in hydrogen-bonded materials [3J El U\ ■ 

Similar potentials arise in the field theory contexts as well with application to domain walls in carbon 
nanotubes and fullerenes |5| . For triangular or hexagonal symmetry crystals two different order parameters 
(e.g. strain and shuffle) can couple with each other being described by an asymmetric double well [HI HUj. 
We consider this situation in detail here and obtain several exact domain wall solutions of a coupled 
asymmetric double well model. Note that in the uncoupled limit, the well known solutions are the (i) kink 
solution (when the condition b\ = Aa\C\ is satisfied, see below) and (ii) the pulse solutions around the 
false vacuum both above and below T c . In this paper we generalize the well known kink solution of the 
uncoupled asymmetric double well model by obtaining a one parameter family of kink solutions of the same 
model. We also obtain a few other solutions in the uncoupled limit, which to the best of our knowledge 
are not known in the literature. 

Lattice dynamical models of entropy driven first order transitions modeled by asymmetric double well 
potentials have been studied previously using both molecular dynamics and mean field theory [111 I12| . 
Exactly solvable asymmetric double well potentials are also known 13... There are many other physical 
contexts in which either the quantum mechanical eigenvalue problem or the domain wall solutions of the 
asymmetric double well potentials have been studied ^1 EGO EE3 EEH EH HO] • Recently, we have obtained 
a large number of exact periodic domain wall solutions for the coupled (ft 4 21j and the coupled (ft 6 [22] 
models. Here we apply the same procedure to coupled asymmetric double well models. 

The paper is organized as follows. In the next section (Sec. II) we present the coupled model and 
obtain the relevant equations of motion. In Sec. Ill we obtain several solutions including the kink lattice 
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("bright-bright" soliton) solution at the transition temperature (T = T c ). Section IV deals with the pulse 
lattice ("dark-dark") solutions both above and below T c . We also obtain the energy of the soliton solutions 
and the asymptotic interaction between the solitons. In Sec. V we consider a different variant of the model 
and obtain a few solutions of the corresponding coupled model. Finally, in Sec. VI we conclude with 
remarks on physical relevance of these solutions and their comparison with the coupled 6 and cf> A models. 



2 Model 

We consider a potential that is asymmetric in the two fields due to a linear-quadratic coupling in addition 
to a biquadratic coupling. The potential we consider is given by (<q > 0, c 2 > 0) 

V = (ai0 2 - b 1( p 3 + ckP 4 ) + (a 2 ip 2 - 6 2 V 3 + c 2 V> 4 ) + #V> 2 + e^V , (1) 

where aq^, 61,2 , ci,2, d, e are material or system dependent parameters. The static field equations that follow 
from here are 

4> xx = 2ai0 - 36i0 2 + 4ci0 3 + d^ 2 + 2e0V> 2 , 

ipxx = 2a 2 ip - 3b 2 ip 2 + 4c 2 V> 3 + 2#V + 2e(p 2 if) • (2) 

Observe that as long as d / 0, the two field equations are asymmetric and hence bright-dark and dark-bright 
solitons would be distinct. 

For the standard uncoupled model (d = e = 0) one usually takes oq > as it then corresponds to a 
model for first order transition. The sign of b\ merely decides if the other minimum is to the right or left 
of <f> = 0. Hence without any loss of generality, one usually restricts to b\ > and we will primarily stick 
to this choice (i.e. aq, 61, c\ > 0) in this paper except in Sec. V. In this case it is easy to show that while 
4> = is the only minimum in case b\ < (32/9)Gqci, for (32/9)aiCi < b\ < 4(qci, ^ = is the absolute 
minimum while 4> = (fr c is the local minimum, whereas for b\ > 4cqci, the opposite is true. Here 



8ci 

while the local maximum of the potential is always at (j) m and is given by 



[36i + ^{Vb\ - 32aici)] 

! (3) 



[3&i - J{%\ -32aici) 



8^ • (4) 
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At b\ = Aa\C\ we have two degenerate minima at (j) = and <f) = <f) c . For the uncoupled model, it is well 
known that while at T = T c (i.e. b\ = 4aiCi), one has a kink solution 012], for T > T c as well as for 
T <T C one has a pulse solution around the local minimum [2j , in a way related to the decay of the false 
vacuum. 

It is perhaps not well appreciated that even for a\ < 0, the uncoupled model corresponds to a model 
for first order transition. In particular, in that case while <j) = is always the local maximum, for b\ > 0, 
(pd is the absolute minimum while (p C 2 is the local minimum, while for b\ < 0, (p C 2 and 4> c \ interchange 
their roles. Here 4> c i^2 are defined by 



[3&i + J(9bj + 32aici)] 
bcl = 8^ ' (5) 



[36i- J(96?+32oici)] 

<Pc2 = ^ • (6) 

At b\ = they are degenerate (the usual 4 — 4> 2 double well potential). 

In this paper, we primarily focus on the case when a\ > and without any loss of generality take b\ > 0. 
Needless to say that from the consideration of stability, c\ is always taken to be positive throughout this 
paper. However, in Sec. V we also write down a few solutions of the coupled model in case a\ < 0. 

We now write down the periodic soliton solutions of the coupled continuum asymmetric 4 models. 
We will also consider the solutions in the (physically interesting) case of 62 = C2 = 0. We will see that, 
corresponding to the uncoupled limit of b\ = Aaic\ one has ten periodic soliton solutions in the coupled 
case which in the hyperbolic limit correspond to three bright-bright, two dark-dark and four dark-bright 
and bright-dark solutions, while one solution becomes a constant in this limit. In the physically interesting 
case of 62 = C2 = 0, only six out of these ten solutions survive. On the other hand, for T > (<) T c , we are 
able to obtain only one solution of this coupled model which, depending on the value of the parameters, 
is related to the pulse solution around the false vacuum (local minimum), both above and below T c . It is 
interesting to note that some of the solutions exist only in the coupled model but not in the uncoupled 
limit. 

For static solutions the energy is given by 



E 



dx , (7) 



where the limits of integration are from — oo to oo in the case of hyperbolic solutions (i.e., single solitons) 
on the full line. On the other hand, in the case of periodic solutions (i.e. soliton lattices), the limits are 
from —K(m) to +K(m) or from —2K(m) to +2K(m) depending on the period of the corresponding elliptic 
function. Here K{m) [and E(m) below] denote the complete elliptic integral of the first (and second) kind 
Using equations of motion, one can show that for all of our solutions 



1 fdd)\ 2 1 /d-if;\ 2 ' 



2 \dx J 2 \dx J 



+ C, (8) 



where the constant C in general varies from solution to solution. Hence the energy E = E — J C dx is 
given by 

2 /«7„/,\ 2" 



E = E — Cdx 



d±Y + fd± 
dx J \dx 



dx. (9) 



Below we will give explicit expressions for energy in the case of several of the periodic solutions (and hence 
the corresponding hyperbolic solutions). In each case we also provide an expression for the constant C. 

3 Solutions Corresponding to b\ = Aa\C\ (in the Analogous Uncoupled 
Case) 

In this and the next section we consider the model with a\,b\,c\ > 0. We look for the most general 
periodic solutions in terms of the Jacobi elliptic functions sn(x,m), cn(x,m) and dn(x,m) [23] which in 
the decoupled limit correspond to b\ = Aa\c\ (i.e. T = T c ). 
Solution I 

It is not difficult to show that 

4> = F + Asn[D(x + x ),m] , ip = G + Bsn[D(x + x ), m] , (10) 

is an exact solution provided the following eight coupled equations are satisfied 

2aiF - 36iF 2 + 4ciF 3 + dG 2 + 2eFG 2 = , (11) 

2a x A - §b x AF + \2c x F 2 A + 2BdG + 2eAG 2 + AeFBG = -(1 + m)AD 2 , (12) 
- 3hA 2 + 12ciF^l 2 + dB 2 + 2eFB 2 + AeABG = , (13) 



2 Cl A 2 + eB 2 = mD 2 , (14) 

2a 2 G - 3b 2 G 2 + 4c 2 G 3 + 2dFG + 2eGF 2 = , (15) 

2a 2 B - Qb 2 BG + \2c 2 G 2 B + 2dAG + 2dFB + AeAFG + 2e£F 2 = -(1 + m)BD 2 , (16) 

- 3b 2 B 2 + 12C2GB 2 + 2dAB + 2eGA 2 + AeABF = , (17) 

2c 2 B 2 + eA 2 = mD 2 . (18) 

Here A and -B denote the amplitudes of the "kink lattice", F and G are constants, D is an inverse 
characteristic length and xq is the (arbitrary) location of the kink. Five of these equations determine 
the five unknowns A, B, D, F, G while the other three equations, give three constraints between the eight 
parameters 01,2, &i,2> d, e. In particular, A and B are given by 

2 _ mD 2 (2c 2 - e) 2 _ mD 2 (2 Cl - e) 

(4c lC2 -e 2 ) ' (4c lC2 -e 2 ) ' [ V) 

It may be noted here that in case both F, G = then no solution exists so long as a\ > 0. However, a 
solution exists in case G = 0, F 7^ or when F = 0, G 7^ 0, both of which we now discuss one by one. 
G=0, F/0: 

In this case A, B are again given by Eq. (|19|) while F and D are given by 



3&i - J{9bj - 32aici)] 



D 2 



(96f - 32aici 



8(1 + m)ci 

Further there are three constraints given by 



8ci 

61 - ^/(96f -32aici) 

1 • (20) 



— 1", 2eF) = 1/ 4 , (6&ici - 4cid - 3e6i)F = 2(2a 2 ci + 4a lCl - eoi) , 

3o 2 y (2c 2 - e) 

27b 2 2 {b 1 -Ac l F)=A(d + 2eF) 2 . (21) 

It is worth noting that the solution exists only when the value of F corresponds to <p m as given by Eq. 0. 
From one of the constraint [Eg. ([21)1] it follows that b\ > A.c\F, i.e. b 2 > Aa\Ci. 
Special case of e 2 = 4cic 2 
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One can show that the solution (jlUj) exists even in case e 2 = Ac\c 2 . It turns out that such a solution 
exists only if 

2ci = 2c 2 = e . (22) 
In this case, A and i? are determined from the relations 

777 /~)^ 

A 2 + B 2 = , 3b 2 B = 2(d + 2eF)A, (23) 

e 

while D, -F are still given by Eq. Q2UJI and the two constraints are 

276| (6i - 2eF) = 4(d + 2eF) 2 , ai + a 2 + dF = . (24) 

Interesting case of b 2 = c 2 = 

Finally, let us discuss the physically interesting case of b 2 = c 2 = 0. In this case, solution (|10[) with 
G = 0,F ^ exists only if e > 2ci. In particular, in this case Eq. (|1U|) with G = is a solution to the 
field Eq. © provided 

2 _ mat „ 2 _ (e - 2ci)mai 2 _ a i 77 _ A^i" 



(l + m)e' (l + m)e 2 ' (1+m) V 4ci V ; 

and further if the following three constraints are satisfied 

6 2 = 4oiCi, d 2 = 2e(oi + 2a 2 ) , M + 2aie = 0. (26) 

F=0, G/0: 

Let us now discuss the solution in case F = but 0. In this case solution (|1U|) exists only if cZ = 
and e < 0. We therefore write e = — |e| and A, B are now given by 

A2 _ mD 2 (2c 2 + \e\) ^ b2 _ mD 2 (2 Cl + \e\) ^ 
(4cic 2 - e 2 ) ' (4cic 2 - e 2 ) 

while G and I? are given by 



[3b 2 - yj(9ti* - 32a 2 c 2 )} 

G = Wm — 



D 2 



(9bl - 32a 2 c 2 ) 



8c 2 

b 2 - i/(96| -32a 2 c 2 ) 

1 • (28) 



8(1 + m)c 2 



Further there are three constraints given by 



3M = -4\e\BG , 3B 2 {b 2 - 4c 2 G) = -2\e\GA 2 , 
3(2c 2 + \e\)b 2 G = 2(2 ai c 2 + 4a 2 c 2 + \e\a 2 ) . 



(29) 



Note that one of the constraint implies that b 2 < 4c 2 G which in turn implies that b\ < 4a 2 c 2 . 

It is easily checked that no solution of the form (|lUj) exists when either e 2 = 4c±c 2 or if 62 = c 2 = 0. 



m=l 



In the limiting case G = 0, F ^ 0, m = 1, we have a bright-bright soliton solution given by 
4> = F + A tanh[D(x + x )] , tp = B tanh[D(x + x )] , 



(30) 



with A, B, F and D given by Eqs. (|19jl and (|2Uj) with m = 1 while the three constraints are again given 
by Eq. (|21|) . A similar solution also exists in case F = 0,G/0. 
Uncoupled Case 

For completeness, it may be worthwhile to write down the solution of the uncoupled field Eq. ©, i.e. 
when d = e = 0. It is easily shown that in this case the solution (|1U|) reduces to 

2m 



^ 1 + 

4ci V (1 + m) 



sn 



«i 



(x + x Q ),m 



(31) 



(1 + m) 

and one is at the transition temperature T = T c since the parameters a\,bi,cx satisfy 6 2 = 4a\Ci. At 
m = l, this reduces to the well known kink solution 



1 + tanh(W^-[x + a?o]) 



(32) 



Energy: Corresponding to the periodic solution (|l(Jj) with G = the energy is same irrespective of whether 
F = 0, G ^ or G = 0, F ^ or if both of them are zero. Only the value of C is different for the different 
cases. For example, for G = 0,F / the energy E and the constant C are given by 

2(A 2 + B 2 )D. 



E 



3m 



-[(1 + m)E(m) - (1 - m)K(m)] , 



1 



C = F 2 [ ai - biF + Cl F 2 } - -(A 2 + B 2 )D 2 . 
On the other hand, in case F = 0, G ^ then C is given by 

C = G 2 [a 2 - b 2 G + c 2 G 2 } - X -{A 2 + B 2 )D 2 . 



(33) 



(34) 



On using the expansion formulas for E(m) and K(m) around m = 1 as given in [23] 



K(m) = In 



+ 



(1 — m) 



E(m) = 1 + 



y/l-mj ' J 
(1 — m) 



In 



In 



m 

1 
2 



+ ..., 



+ ... 



(35) 



(36) 



— VTLj 

for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (bright-bright) soliton solution [Eq. Q30J1 ] plus the interaction energy. We find 



E = Ekink + E int = (A 2 + B 2 )D 



4 (1-m) 

- + - 

3 3 



The interaction energy vanishes at exactly m = 1, as it should. 
Solution II 

A different type of solution (i.e. a pulse lattice) given by 

4> = F + Acn[D(x + xo), m] , ip = G + Bcd[D{x + xo),m] , 

is an exact solution provided the following eight coupled equations are satisfied 

2aiF - 36iF 2 + 4ciF 3 + dG 2 + 2eFG 2 = , 

2 ai A - GbiAF + 12 Cl F 2 A + 2BdG + 2eAG 2 + AeFBG = (2m - l)AD 2 , 
- 3hA 2 + UciFA 2 + dB 2 + 2eFB 2 + AeABG = , 
2c\A 2 + eB 2 = -mD 2 , 
2a 2 G - 3b 2 G 2 + 4c 2 G 3 + 2dFG + 2eGF 2 = , 
2a 2 B - Qb 2 BG + \2c 2 G 2 B + 2dAG + 2dFB + AeAFG + 2eBF 2 = (2m - l)BD 2 , 
- 3b 2 B 2 + \2c 2 GB 2 + 2dAB + 2eGA 2 + AeABF = , 
2c 2 B 2 + eA 2 = -mD 2 . 



(37) 



(38) 



(39) 

(40) 
(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



Notice that two of these equations are meaningful only if e < since from stability considerations, c\ > 
0, c 2 > 0. Hence, we write e = — lei. 



Five of these equations determine the five unknowns A, B, D, F, G while the other three equations, give 
three constraints between the eight parameters 0,^2, 61,2, ci,2, d, e. In particular, A and B are given by 

2 = mD 2 (2c 2 + \e\) 2 mD 2 (2 Cl + \e\) 

( e 2-4 Cl c 2 ) ' (e 2 -4c lC2 ) • 1 j 

Unlike the previous case, it turns out that in this case a solution exists both when F = G = 0, G = 
or when F 7^ or if F = 0, G 7^ and we discuss all three cases one by one. 
F=G=0 

In this case A, B are again given by Eq. (|47|) while D, the characteristic inverse length, is given by 

(2m — 1) 

while the three constraints are 



Since ai > 0, such a solution exists only if m > 1/2. It is easily checked that such a solution (with 
F = G = 0) does not exist in case either 62 = ci = or if e 2 = 4ciC2- 

What happens if m = 1/2? It is easily checked that in that case (|38j) with F = G = is still a solution 
provided ai = a 2 = 0, Z) is undetermined, A, B are given by Eq. (|47jl and the two remaining constraints 
are given by Eq. (|49|). 
G=0, F/0: 

In this case 1)38(1 is a solution with .A, B again given by Eq. (|47|) while i 7 and I? are given by 



[3&i + t/(96? - 32oici)] 

F = 4>c = 

8ci 



, 96 2 - 32aici 3&i + J (9b 2 - 32aici 

Z? = - ± . (50) 

8(2m-l)ci v ; 

Further there are three constraints given by 

(2|e|F-d) (2 Cl + |e|) a ^\ \h \tp ofo j.a xl I ^ 

^771 — F TT = 7^ TTTT ' ( 6& i c i ~ 4c i« + 3 e oJF = 2(2a 2 ci +4aiCi + e ai , 

3(4ciF - bi) (2c 2 + |e|) 

276|(4ciF-fei) = A(2\e\F - d) 2 . (51) 
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What happens if m = 1/2? It is easily checked that in that case (|38|) with G = is still a solution 
provided A, B are given by Eq. (|47|) . -D is undetermined while the constraints are 

3&iF = 4ai , 96? = 32a lC i , 3b 2 B = 2dA , (2\e\F - d)B 2 = Ac x FA 2 , a 2 + dF = \e\F 2 . (52) 

Interestingly enough, there is also a solution in case m < 1/2. In particular, it is easily checked that 
(|38|) is a solution with A, B again given by Eq. ()47|) while F and D are now given by 



[36i-^/(96f-32aici)] 

F = 4>rn = ~ 

8ci 



., _ v '(96? - 32a lCl ) [36i - ^/(96? - 32a lCl )] 
D " 8(1 - 2m) Cl • (53) 

Further the three constraints are again given by Eq. (|51j). 
Interesting case of b 2 = c 2 = 

Finally, let us discuss the physically interesting case of b 2 = c 2 = 0. It is easy to see that now a solution 
of the form (|38[) with G = is possible only if m < 1/2. In particular, such a solution exists if 

2 _ mD 2 2 _ m£> 2 (2ci + |ej) 2 _ gl 
^ ~TT' " i 2 ^ ' " (1 - 2m) ' (M) 

and if further the following constraints are satisfied 

6 2 = 4aici, dbi = 2\e\a 1 , d 2 = -2|e|(oi + 2a 2 ) . (55) 

Thus with ai > 0, a solution of the form (|38|1 with b 2 = c 2 = can exist only if m < 1/2 and a 2 < 0. 

It is easily checked that no solution is however possible in the special case of Ac\c 2 = e 2 . 
F=0, G/0: 

In this case (|38|) is a solution only if d = with A, B again given by Eq. (|47|) while G and are given 

by 



[36 2 + J(9^ - 32a 2 c 2 )} 

G = xp c = 

8c 2 



' 96 2 - 32a 2 c 2 36 2 + J (9b 2 - 32a 2 c 2 )} 

D 2 = ^ . (56) 

8 2m - 1 c 2 
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Further there are three constraints given by 

36iA = -4\e\BG, {Ac 2 G - b 2 )B 2 = 2\e\GA 2 , 2oi - 2\e\G 2 = (2m - l)D 2 . (57) 

Note that this solution exists only if m > 1/2. 

What happens if m = 1/2? It is easily checked that in that case (|38|) with F = is still a solution 
provided j4, 5 are given by Eq. (|47|). -D is undetermined and the constraints are 

36iA = -4|e|BG , a x = \e\G 2 , b 2 B 2 = A\e\GA 2 , 36 2 G = 4a 2 , 9^ = 32a 2 c 2 . (58) 

Interestingly enough, there is also a solution in case m < 1/2. In particular, it is easily checked that 
()38|) is a solution with A, B again given by Eq. (|47jl while G and D are now given by 



[36 2 - J(96| - 32a 2 c 2 )] 

G = i' m = 1 5 

8c 2 



^(961 - 32a 2 c 2 ) [3fe 2 - ^/(9j - 32a 2 c 2 )] 
8(l-2m)c2 ' 

The constraints are again given by Eq. (|57j). Note that one of the constraint can only be satisfied if 

b\ < 4a 2 c 2 . Even in this case, no solution of the form (|3*%|) with F = 0,G/0 exists either when e 2 = 4cic 2 

or when 6 2 = c 2 = 0. 

m=l 

In this limiting case we have a dark-dark soliton solution given by 

4> = F + Asech[D(x + x )} , 4) = G + Bsech[D(x + x )} , (60) 

where F and/or G may or may not be zero depending on the solution that we are considering. 

It is worth emphasizing that while solution (|38[) exists in the coupled case, there is no analogous solution 
to the corresponding uncoupled model, i.e. no solution is possible in the case of the uncoupled field Eq. 
(J2J) i.e. when d = e = 0, so long as c\ > 0. 

Energy: Corresponding to the periodic solution (]38[) the energy is same but only the value of C changes 

depending on if F = G = or G = 0, F / or F = 0, G / 0. For example, if G = 0, F ^ then the 

energy E and the constant C are given by 

2( A 2 + B 2 )D 

E = -i '— [(2m - l)E(m) + (1 - m)iC(m)] , 

3m 
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C = F 2 [ ai - biF + ciF 2 ] - 1(1 - m)(A 2 + B 2 )D 2 . (61) 



On the other hand, if F = 0, G ^ then C is given by 



C = G 2 [a 2 - b 2 G + c 2 G 2 } - hi - m)(A 2 + B 2 )D 2 , (62) 

while if F = G = then C is simply given by 

C = -^(1 - m)(A 2 + B 2 )D 2 . (63) 

On using the expansion formulas for E{m) and K{m) around m = 1 as given by Eqs. (|35|) and ()36|) . 
for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-dark) soliton solution [Eq. (|60|) ] plus the interaction energy. We find 



E = Ep U i se + Ei nt = (A 2 + B 2 )D 



2 5(1 -m) M . 4 

; + (1 - m)ln(- 



(64) 



.3 6 y/1 — m' 

Note that this solution exists only when e < and e 2 > ^c\c 2 - The interaction energy vanishes at exactly 
m = 1, as it should. 
Solution III 

Yet another pulse lattice solution given by 

cf) = F + Adn[D(x + x ),m], tp = G + Bdn[D(x + x ), m] , (65) 

is an exact solution provided the following eight coupled equations are satisfied 

2 ai F - 3&iF 2 + 4ciF 3 + dG 2 + 2eFG 2 = , (66) 

2 ai A - 6&i AF + Ylc x F 2 A + 2BdG + 2eAG 2 + AeFBG = (2 - m)AD 2 , (67) 

- 36iA 2 + 12ciF^ 2 + dB 2 + 2eFB 2 + 4eABG = , (68) 

2ci^ 2 + eB 2 = —D 2 , (69) 

2a 2 G - 3b 2 G 2 + 4c 2 G 3 + 2dFG + 2eGF 2 = , (70) 

2a 2 B - 6b 2 BG + l2c 2 G 2 B + 2dAG + 2dFB + 46^(7 + 2eSF 2 = (2 - m)BD 2 , (71) 
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- 3b 2 B 2 + 12c 2 GB 2 + 2dAB + 2eGA 2 + AeABF = , (72) 

2c 2 B 2 + eA 2 = —D 2 . (73) 

Notice that two of these equations are meaningful only if e < 0, since from stability considerations, 

c\ > 0, c 2 > 0. Hence we write e = — |e|. 

Five of these equations determine the five unknowns A, B, D, F, G while the other three equations, give 

three constraints between the eight parameters ai t 2, b\ 2, c\ 2, d, e - I n particular, A and B are given by 

2 = D 2 (2c 2 + \e\) 2 D 2 (2c 1 + |e|) 

( e 2-4 Cl c 2 ) ' (e 2 -4c lC2 ) • {<V 

Akin to the solution II, in this case too, a solution exists when i ? = G = 0orG = 0,F^0or when 

F = 0, G 7^ and we discuss all the cases one by one. 

F=G=0 

In this case A, B are again given by Eq. (|74|) while D is given by 

D 2 = , (75) 

(2 - m) v ; 

and the three constraints are 



It may be noted that such a solution with F = G = will not exist in case either b 2 = c 2 = or if 
e 2 = 4cic 2 . 
G=0, F/0: 

In this case A, B are again given by Eq. (|74|) while F and D are given by 



3&i + V(96f - 32oici) 
F = 6 r = 



8ci 



, 96? - 32aici 3&i + J (9b 2 - 32a lCl )] 

D = ^ ; y . (77) 

8(2 - m)ci v ' 

Further there are three constraints given by 

<2\e\F — d) (2a + lei) 

q / "7 ' , \ = , , (66ici - Add + 3|e|6i)F = 2(2a 2 ci + Aa lCl + |e|oi) , 
3(4ciF - Oi) (2c 2 + |e|) 

276^(4ciF-6i) = A(2\e\F - d) 2 . (78) 
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It is worth noting that since b\ — 4ciF < 0, hence F > d/2\e\. It is easily checked that no solution is 

however possible in case b 2 = c 2 = or if 4ciC2 = e 2 . 

m=l 

In this limiting case we have the same dark-dark soliton solution as given by Eq. (|6Uj) . 
F=0, G/0: 

Solution (|65|) with F = 0, G / is only possible if d = 0. In this case A, B are again given by Eq. (|74j) 
while G and D are given by 



36 2 + V (96 2 _ 32a2C2 ) 
G = ip c = 



8c 2 



;(962 - 32a 2 c 2 ) [36 2 + Vgg - 32a2C 2 )] 
8(2 -m)c 2 ' 

Further there are three constraints given by 

3hA = 4\e\BG, 2 ai - 2\e\G 2 = (2 - m) D 2 , 3(4c 2 G - b 2 )B 2 = 2\e\GA 2 . (80) 

Again, it is easily checked that no solution is however possible in case 62 = C2 = or if 4ciC2 = e 2 . 

It may be noted that, as in the previous case, while there exists a solution (|65j) in the coupled case, no 
such solution exists in the corresponding uncoupled case, i.e. no solution exists in the case of the uncoupled 
field Eq. (J2J) i.e. when d = e = 0, so long as c\ > 0. 

Energy: Corresponding to the periodic solution (|65|). the energy is same but only C is different depending 
on if F = G = 0, or G = 0, F / or F = 0, G / 0. For example, with G = 0, F ^ 0, the energy E and the 
constant C are given by 

E = 2 ( A2 + - B2 ) g [(2 - m)E(m) - (1 - m)K{m)} , 
3 

C = F 2 [ ai - b x F + Cl F 2 ] + 1(1 - m )(A 2 + S 2 )!) 2 . (81) 
On the other hand, if -F = 0, G 7^ 0, then C is given by 

C = G 2 [a 2 - b 2 G + c 2 G 2 } + 1(1 - m)(A 2 + S 2 )!) 2 > (82 ) 
while if F = G = then C is simply given by 

C = ^(1 - m)(A 2 + B 2 )D 2 , (83) 
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On using the expansion formulas for E(m) and K(m) around m = 1 as given by Eqs. (|35|) and ()36j) . 
for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-dark) soliton solution [Eq. (|6U|) ] plus the interaction energy. We find 



E = E pu i se + E int = (A 2 + B 2 )D 



2 (l-m) ( 4 

(1 — m) in 



(84) 



_3 2 v ' \^r~m, 

Note that this solution exists only when e < 0,e 2 > 4cic 2 . The interaction energy vanishes at exactly 
m = 1, as it should. Unlike the dark-dark cn — cn and dn — dn solutions, it turns out that the (mixed) 
dark-dark periodic soliton solutions of the form dn — cn and cn — dn, do not exist. 
Solution IV 

We now discuss two bright-dark and two dark-bright periodic soliton solutions. We will see that for 
these four solutions G is necessarily zero while F is necessarily nonzero. Let us discuss these solutions one 
by one. In particular, it is easily shown that 

4> = F + Asn[D(x + Xo),m] , V = G + Bcn[D{x + x ), m] , (85) 

is an exact solution provided 

G = , b 2 = , d + 2eF = , (86) 
and the following six coupled equations are satisfied 

2aiF- 36iF 2 + 4ciF 3 = 0, (87) 

2oi - 661F + 12ciF 2 + 2eB 2 = -(1 + m)D 2 , (88) 

- 36iA 2 + \2ciFA 2 = 0, (89) 

2ci^ 2 - eB 2 = mD 2 , (90) 

2a 2 + 2dF + 2e(A 2 + F 2 ) = (2m - 1)£> 2 , (91) 

2c 2 B 2 - eA 2 = -mD 2 . (92) 

Three of these equations determine the three unknowns A,B,D while the other three equations, give three 

constraints between the seven parameters ai^, 6i, ci^, d, e. In particular, A and B are given by 

A2 _ mD 2 (e - 2c 2 ) ^ ^ 2 _ mZj 2 (2 Cl -e) ^ 
(e 2 - 4cic 2 ) ' (e 2 - 4cic 2 ) 
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while the inverse characteristic length D is given by 

D 2 = (e 2 -4 ClC2 ) ai 

[(1 + m)(e 2 -4cic 2 ) + 2em(2ci - e)] ' 

Further, the three constraints are 

6 2 = 4aici, 2aie + d6i = 0, (2m - 1)D 2 = 2a 2 + 2e^ 2 + dF . (95) 

m=l 

In this limiting case we have a bright-dark soliton solution given by 

4> = F + A tanh[I>(x + x )] , tp = Bsech[D(x + x )] , (96) 

with A, B and D given by Eqs. (|93ft and (|94jl with m = 1 while the three constraints are again given by 
Eq. 

Special case of e 2 = 4ciC2 

One can show that the solution (jHSJ) exists even in case e 2 = 4cic 2 . It turns out such a solution exists 
only if Eqs. (|22|) and (|86|) are satisfied. In this case, A, B and D are given by 

A 2 = t mfl i ~ 2 (! ~ ra)a 2 ] ^2 = ~ 2(1 + m)q 2 ] D 2 = ^2 
2em 2em m ' 

while the constraints are 

d = -h, b\ = 4oiCi . (98) 

Interesting case of 6 2 = c 2 = 

Finally, let us discuss the physically interesting case of 6 2 = c 2 = 0. In this case Eq. (|85|) with G = 
is a solution to the field Eqs. @ provided 

a2 = m{d 2 - 4ea 2 ) ^ 2 = (2ci - e)(d 2 - 4ea 2 ) ^ 2 = {d 2 - 4eo 2 ) = _d_ 

2e 2 ' 2e 3 ' 2e 2e ' 1 j 

and further if the following three constraints are satisfied 

6^ = 4aici , (d 2 - 4ea 2 )[4mci + (1 - m)e] = 2aie 2 , M + 2aie = 0. (100) 
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Energy: Corresponding to the periodic solution (|85[) the energy E and the constant C are given by 

E = — (\(2 - m)mB 2 + (1 + m)A 2 ]E{m) - (1 - m)(A 2 + 2B 2 )K(m)) , 
3m V / 

C = F 2 [ ai - biF + Cl F 2 } - ^A 2 D 2 + B 2 [a 2 + dF + eF 2 + c 2 B 2 } . (101) 

On using the expansion formulas for E(m) and K(m) around m = 1 as given by Eqs. (|35j) and Q36JI . 
for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (bright-dark) soliton solution [Eq. I)96j) ] plus the interaction energy. We find 



F — E so n ton + Ei n t — D 



-(2A 2 + B 2 ) + (1 m) (2A 2 + 3B 2 ) - (1 - m)B 2 In ( ' 



(102) 



.3 6 \ \J\ — m, 

Note that this solution exists when either (i) 2c\ > e > 2c 2 and e 2 > kc\c 2 or if (ii) 2c 2 > e > 2ci and 
Ac\c 2 > e 2 or if (iii) c 2 = 0, 2ci > e. The interaction energy vanishes at exactly m = 1, as it should. 
Solution V 

It is easy to show that another such (kink- and pulse-like) solution is 

(f) = F + Asn[D(x + x ),m] , ip = G + Bdn[D(x + x ),m] . (103) 
This is an exact solution provided 

G = 0, b 2 = 0, d + 2eF = 0, (104) 

and the following six coupled equations are satisfied 

2a x F- 3hF 2 + 4ciF 3 = 0, (105) 

2ai - 6hF + 12ciF 2 + 2eB 2 = - (1 + m)D 2 , (106) 

- 36iA 2 + 12ciFy4 2 = , (107) 

2ci^ 2 - eraB 2 = mD 2 , (108) 
A 2 

2a 2 + 2dF + 2e( — + F 2 ) = (2 - m)D 2 , (109) 
m 

2mc 2 B 2 - eA 2 = -mD 2 . (110) 



18 



Three of these equations determine the three unknowns A, B, D while the other three equations, give three 
constraints between the seven parameters ai t 2, &i, c± t 2, d, e. In particular, A and B are given by 

2 _ mD 2 (e-2c 2 ) 2= D 2 {2 Cl -e) 

( e 2-4c lC2 ) ' a ( e 2_ 4ciC2 )' V 11 ) 



while D is given by 



& = (e 2 -4 Cl c 2 )ai 

[(l+m)(e 2 -4cic 2 ) + 2e(2ci - e)] ' V ; 



Further, the three constraints are 

b? = 4oic ls 2aie + d6i=0, (2 - m)mD 2 = 2a 2 m + 2eA 2 + mdF . (113) 

m=l 

In this limiting case we have a bright-dark soliton solution given by Eq. (|96|) with A, -B and D given 
by Eqs. Q93|) and 1)94(1 with m = 1 while the three constraints are again given by Eq. (|95|) . 
Special case of e 2 = 4cic 2 

One can show that the solution (|1U3|) exists even in case e 2 = 4cic 2 . It turns out such a solution exists 
only if Eqs. (|22|) and (|1U4)I are satisfied. In this case, A, B and -D are given by 

2 = m[oi + 2(1 -m)q 2 ] 2 = [gi - 2(1 + m)a 2 ] 2 = 

2e 2e 2 ' 1 j 

while the constraints are 

d = -6i , 6? = 4oici • (115) 

Interesting case of b 2 = c 2 = 

Finally, let us discuss the physically interesting case of b 2 = c 2 = 0. In this case Eq. (|1U3|) with G = 
is a solution to the field Eq. (J2J) provided 

A 2 = (d 2 J -4eo 2 ) g2 _ (2ci-e)(d 2 -4ea 2 ) j2 _ ((i 2 - 4eo 2 ) p = _± (nQ) 
2e 2 ' 2me 3 ' 2me ' 2e ' 1 ' 

and further if the following three constraints are satisfied 

b\ = 4oici , (d 2 - 4ea 2 )[4ci - (1 - m)e] = 2maie 2 , 6id + 2aie = 0. (117) 
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Energy: Corresponding to the periodic solution ()103l) the energy E and the constant C are given by 

E = — ([{2m - l)B 2 + (1 + m)A 2 ]E(m) - (1 - m)(A 2 - B 2 )K(m)\ , 
3m \ J 

C = F 2 [ ai - biF + ciF 2 } - ^A 2 D 2 + B 2 [a 2 + dF + eF 2 + c 2 B 4 } . (118) 

On using the expansion formulas for E(m) and K(m) around m = 1 as given by Eqs. (|35j) and ()36|) . 
for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (bright-dark) soliton solution [Eq. I)96j) ] plus the interaction energy. We find 



E — E so n ton + Ei n t — D 



-(2A 2 + B 2 ) + (1 m) (2A 2 - 5B 2 ) + (1 - m)B 2 In ( ' 



(119) 



.3 6 \ -y/1 — m, 

Note that this solution exists when either (i) 2c\ > e > 2c 2 and e 2 > kc\c 2 or if (ii) 2c 2 > e > 2ci and 
4ciC2 > e 2 or if (iii) c 2 = 0, 2c\ > e. The interaction energy vanishes at exactly m = 1, as it should. 
Solution VI 

We now discuss two periodic solutions which at m = 1 reduce to dark-bright soliton solutions. In 
particular, it is easily shown that 

(f) = F + Acn[D(x + x ),m] , if) = G + sn[D(x + x ),m] , (120) 

is an exact solution provided 

G = 0, b 2 = 0, d + 2eF = 0, (121) 
and the following six coupled equations are satisfied 

2a x F - 3hF 2 + 4ciF 3 = 0, (122) 

2ai - 6hF + 12ciF 2 + 2eB 2 = (2m - l)D 2 , (123) 

- 3biA 2 + 12ciF,4 2 = , (124) 

2ci^ 2 - eB 2 = -mD 2 , (125) 

2a 2 + 2dF + 2e(^ 2 + F 2 ) = -(1 + m)D 2 , (126) 

2c 2 B 2 - eA 2 = mD 2 . (127) 
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Three of these equations determine the three unknowns A, B, D while the other three equations, give three 

constraints between the seven parameters a% t 2, b\, c\ t 2, d, e. In particular, A and B are given by 

2 _ mD 2 (e-2c 2 ) p2 _ mfl 2 (2 Cl -e) 
A ~ (4 Cl c 2 -e 2 ) ' a - (4c lC2 - e 2) ' {LZti) 



while D is given by 



D 2 = (4c lC2 -e 2 ) Ql 

[(4cic 2 - e 2 ) + 4mci(e - 2c 2 )] 



Further, the three constraints are 

6 2 = 4aici, 2aie + d6i = 0, - (1 + m)D 2 = 2a 2 + 2eA 2 + dF . (130) 

m=l 

In this limiting case we have a dark-bright soliton solution given by 

4> = F + Asedi[D(x + x )], V = B tanh[D(x + x )} , (131) 
with A, i? and D given by Eqs. (|128|) and (j!29j) with m = 1 while the three constraints are given by Eq. 

<nm 

Special case of e 2 = 4cic 2 

One can show that the solution (|120j) exists even in case e 2 = 4cic 2 . It turns out such a solution exists 
only if Eqs. ()22l) and (|121j) are satisfied. In this case, A, B and D are given by 

2 (mai + 2a 2 ) 2 [max - 2(2m - l)a 2 ] 2 2a 2 

^ = o ' B = n ' U = ' ( 162 ) 

2em 2em m 

while the constraints are 

d=-h, b\ = 4aici . (133) 

Finally, let us discuss the physically interesting case of 6 2 = c 2 = 0. In this case, a solution does 
not exist since whereas A 2 > requires that e < but that in turn makes B 2 < 0, since from stability 
considerations, c\ > 0. 

Energy: Corresponding to the periodic solution (|12UI) the energy E and the constant C are given by 

E = — ( [(2m - l)A 2 + (1 + m)B 2 ]E(m) - (1 - m)( J B 2 - A 2 )K(m)] , 
3m V / 

C = F 2 [ ai - hF + ciF 2 } + -(1 - m)^ 2 !) 2 + S 2 [a 2 + + eF 2 + c 2 B 2 } . (134) 
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On using the expansion formulas for E{m) and K(m) around m = 1 as given by Eqs. (|35j) and Q36JI . 
for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-bright) soliton solution [Eq. (|131j) ] plus the interaction energy. We find 



E — EsoUton + Ei n t — D 



l(2B 2 + A 2 ) + (2B 2 - 5 A 2 ) + (1 - m)A 2 In ( ' 



(135) 



.3 ' 6 ' \^J\ — m, 

Note that this solution exists when either (i) 2c\ > e > 2c 2 and 4ciC2 > e 2 or if (ii) 2c2 > e > 2c\ and 
e 2 > 4ciC2- The interaction energy vanishes at exactly m = 1, as it should. 
Solution VII 

Another periodic solution which at m = 1 reduces to the dark-bright soliton solution is given by 

(f) = F + Adn[D(x + x ),m], ip = G + Bsn[D(x + x ), m] , (136) 

provided 

G = 0, b 2 = 0, d + 2eF = 0, (137) 
and the following six coupled equations are satisfied 

2a x F - 3hF 2 + 4ciF 3 = 0, (138) 
2eB 2 

2ai - 661F + 12 Cl F 2 + = (2 - m)D 2 , (139) 

m 

- 36!^1 2 + \2c x FA 2 = , (140) 

2m Cl A 2 - eB 2 = -mD 2 , (141) 

2a 2 + 2dF + 2e(A 2 + F 2 ) = -(1 + m)D 2 , (142) 

2c 2 B 2 - meA 2 = mD 2 . (143) 

Three of these equations determine the three unknowns A, B, D while the other three equations, give three 
constraints between the seven parameters ai2,bi,ci2,d,e. In particular, A and B are given by 

JC-^M, B 2 = "f 2 ' 2 °' - '» , (144) 
(4cic 2 - e 2 ) (4cic 2 - e 2 ) 



while -D is given by 



D 2 = (4cic 2 - e 2 ) ai 

[m(4cic 2 - e 2 ) + 4c a (e - 2c 2 )] V ' 
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Further, the three constraints are 

&?=4o i c 1 , 2 ai e + dbi = 0, - (1 + m)D 2 = 2a 2 + 2eA 2 + dF . (146) 

m=l 

In this limiting case we have a dark-bright soliton solution given by Eq. (|131j) . 
Special case of e 2 = 4ciC2 

One can show that the solution (|136j) exists even in case e 2 = ^c\c 2 . It turns out such a solution exists 
only if Eqs. (122() and (|137j) are satisfied. In this case, A, £> and -D are given by 

2= (a 1 + 2ma 2 ) ^ = [ax - 2(2 - m)a 2 ] 2 = 

2e 2e v ; 

while the constraints are 

d = -h, b\ = 4oiCi . (148) 

Finally, let us discuss the physically interesting case of b 2 = c 2 = 0. As in the previous case, a solution 
does not exist since whereas ^4 2 > requires that e < but that in turn makes B 2 < 0, since from stability 
considerations, c\ > 0. 

Energy: Corresponding to the periodic solution (|136|) the energy E and the constant C are given by 

E = — ([(2 - m)mA 2 + (1 + m)B 2 ]E(m) - (1 - m)(B 2 + 2A 2 )K(m) \ , 
3m V ' J 

1 B 2 

C = F 2 \ai - biF + ciF 2 ] + -(1 - m)A 2 D 2 + -^[ma 2 + mdF + meF 2 + c 2 B 4 ] . (149) 

2 m z 

On using the expansion formulas for E{m) and K{m) around m = 1 as given by Eqs. (|35|) and (|36|) . 
for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-bright) soliton solution [Eq. H131|) ] plus the interaction energy. We find 



E — Egoiitan + Ei n t — D 



1(2B 2 + A 2 ) + (2B 2 + 3A 2 ) - (1 - m)A 2 In ( ' 



(150) 



.3 6 \a/1 — m, 

Note that this solution exists when either (i) 2c\ > e > 2c2 and 4ciC2 > e 2 or if (ii) 2c2 > e > 2ci and 
e 2 > 4ciC2- The interaction energy vanishes at exactly m = 1, as it should. 
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Solution VIII 

We now present three solutions, which to the best of our knowledge were not known before even in the 
uncoupled limit. Two of these solutions, in the uncoupled limit are kink-type while one is a pulse type 
solution. Let us discuss these solutions one by one. 

One of these periodic solutions is given by 

,_ F , Asn[D(x + x ),m] Hsn[D(x + x ), m] 

9 l + Bdn[D{x + x ),m] ' V 1 + Bdn[D(x + x ), m] ' 1 ' 

provided 

B = l, (152) 
and the following eight coupled equations are satisfied 

2aiF - 3biF 2 + A Cl F 3 + dG 2 + 2eFG 2 = , (153) 

2a x A - Q^FA + 12 Cl F 2 A + 2dGH + AeFGH + 2eAG 2 = - A ^ ~™^° 2 , (154) 

- 3hA 2 + \2c x FA 2 + d# 2 + 2eF# 2 + AeAGH = , (155) 

8ciA 2 + 4e# 2 = m 2 D 2 , (156) 

2a 2 G - 3b 2 G 2 + 4c 2 G 3 + 2dFG + 2eGF 2 = , (157) 

2a 2 H - 6b 2 GH + 12c 2 G 2 H + 2dFH + 2dAG + AeFGA + 2etfF 2 = _ H ( 2 -™) D ? ( 158 ) 

- 36 2 # 2 + l2c 2 GH 2 + 2tL4# + 2eCM 2 + 4eAF# = , (159) 

8c 2 # 2 + 4e.4 2 = m 2 L> 2 . (160) 

Five of these equations determine the five unknowns A, H, D, F, G while the other three equations, give 
three constraints between the eight parameters a± t 2, &i,2> c\,2, d, e. In particular, A and H are given by 



4(4cic 2 -e 2 ) ' 4(4cic 2 -e 2 ) 



It turns out that there is no solution to these field equations with F = G = so long as a\ > 0. 
However, there are solutions when either G = 0, F/0orF = 0,G/0 which we discuss one by one. 
G=0, F / 
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In this case A, H are again given by Eq. ()161|) . B = 1 while F and D are given by 



3&i - ^{%\ - 32aici)] 
8c[ 



2 v '(% 2 - 32aici)[36i - (9b( - 32a lCl )] 



D = ^ J . (162) 

4(2 - m)ci V ; 



Further there are three constraints given by Eq. (|21|) . 
Special case of e 2 = 4ciC2 

One can show that the solution Q151j) with G = 0, a\ = &i_F exists even in case e 2 = 4ciC2 provided if 
Eq. (|22|) is satisfied. In this case, A and H can be determined from the relations 

A 2 + H 2 = ^-^-, 3b 2 H = 2(d + 2eF)A, (163) 
4e 

F, D are still given by ()162|) . while the two constraints are given by Eq. (|2"i)) . 
Interesting case of b 2 = c 2 = 

Finally, let us discuss the physically interesting case of b 2 = c 2 = 0. In this case Eq. (|151|) with G = 
is a solution to the field Eq. © provided 

A 2 = f ai , , ff 2 = (e -J Cl)mV , D 2 = t 2a \, F= x f^, (164) 
2e(2-m)' 2e 2 (2-m) ' (2 - m) ' V 4c x ' v ; 

while the three constraints are given by Eq. (|26jl . 

Thus the solution with G = 0,F/0 exists when 2ci > e, 2c2 > e and 4ciC2 > e 2 . 
F=0, G / 

Solution (|15 If) is a solution with F = 0,G / provided d = 0, B = 1 and e < 0, we therefore write 
e = — \e\. In this case A, H are given by Eq. (|161|) while G and D are given by 



[36 2 - V(96g - 32a 2 c 2 ) 

G = Vm = 



8c 2 



, ,/(96 2 - 32a 2 c 2 )[36 2 - J (9b 2 - 32a 2 c 2 )] 

£) = 77 7^ • (165) 

4(2 - m)c 2 y ' 

Further there are three constraints given by 

3hA = -4\e\GH , 3H 2 (b 2 - Ac 2 G) = -2\e\GA 2 , 3(2c 2 + \e\)b 2 G = 2(2 ai c 2 + 4a 2 c 2 + \e\a 2 ) . (166) 
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For F = 0, G 7^ 0, no solution of the form ()151j) exists in case either e 2 = 4ciC2 or if 62 = C2 = 0. 
m=l 

In this limiting case we have a bright-bright soliton solution given by 



F + A tanh 



ai < , ^ 



, ip = i^tanh 



(167) 



with A and given by Eq. ()161|) . D 2 = 2a\,B = 1 while the three constraints are again given by Eq. 
(f2Tj) . Uncoupled Case 

For completeness, it may be worthwhile to write down the solution of the uncoupled field Eq. ©, i.e. 
when d = e = 0. It is easily shown that in this case the solution ()151|) reduces to 

G/£(s + so),m) 



1 + 



m 
2 — m 



sn 1 



(168) 



_ 1 + dn (\/S^( a; + 2; o),™)_ 

and one is at the transition temperature T = T c since the parameters oi,6i,ci satisfy b\ = Aa\C\. At 
m = 1, this reduces to the well known kink solution (|32[). 

Energy: Corresponding to the periodic solution <|151|) . while the energy is same but the value of C is 
different for solutions with either G = 0, F 7^ or F = 0, G 7^ 0. For example, in case G = 0, F 7^ the 
energy £7 and the constant C are given by 



J-2if [1 + dn(y, m)f 
C = [ai - 61F + ciF 2 ]F 2 - i(^ 2 + H 2 )D 2 , 



where y = -D(x + xq). On the other hand, when F = 0, G 7^ 0, the constant C is given by 



(169) 



C = [a 2 - b 2 G + c 2 G 2 ]G 2 - \ {A 2 + # 2 )D 2 . 



(170) 



Unfortunately, we are unable to solve the integral (|169j) analytically and hence are unable to calculate 
the corresponding soliton interaction energy. However, the energy of the corresponding hyperbolic (bright- 
bright) soliton solution [Eq. pbTj) ] is easily calculated. We find 



E 



soliton 



16a 



1 t A2 



18 



[A 1 + H Z )D . 



(171) 
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Solution IX 

We now present a periodic solution, which at m = 1 reduces to a constant but is nontrivial for arbitrary 
m. This periodic solution is given by 

^ _ F | Acn[D(x + x ),m] _ Q | Hcn[D(x + a ), m] 2 

\/l — m + i?dn[Z?(x + xo), m] \/l — m + Bdn[D(x + xo),m] 

provided Eq. (j!52[) holds (i.e. B = 1) and the following eight coupled equations are satisfied 

2aiF - 3&iF 2 + 4ciF 3 + dG 2 + 2eFG 2 = , (173) 

A(2 — m)D 2 

2a x A - QbiFA + 12 Cl F 2 A + 2dGH + AeFGH + 2eAG 2 = i — ^ — , (174) 

- 3hA 2 + \2c x FA 2 + dH 2 + 2eFH 2 + AeAGH = , (175) 

SciA 2 + 4e^ 2 = m 2 D 2 , (176) 

2a 2 G - 3b 2 G 2 + 4c 2 G 3 + 2dFG + 2eGF 2 = , (177) 

2a 2 H - 6b 2 GH + 12c 2 G 2 H + 2dFH + 2dvlG + 4eFG^l + 2eHF 2 = i — — ^ , (178) 

- 36 2 ^ 2 + l2c 2 GH 2 + 2dA/7 + 2eGA 2 + AeAFH = , (179) 

8c 2 H 2 + 4e^ 2 = m 2 L> 2 . (180) 

It is rather remarkable that the eight field equations that we obtain here are identical to those obtained 
for the solution VIII. As a result, all the solutions obtained in that case continue to be true even in this 
case. In particular, as in the previous case, there does not exist a solution to Eq. ()172|) in case F = G = 0. 
On the other hand, when G = 0, F ^ or when F = 0, G ^ 0, the solutions exist with the parameters as 
given for solution VIII. 
m=l 

In this limit, (unlike the solution VIII) the solution (|172[) with G = 0, F ^ goes over to a constant 
solution, given by 

cf) = F + A, ip = H, (181) 

where A, H are given by Eq. ()161j) while F = <j) m with <f> m given by Eq. Similar conclusions are also 
obviously valid for the solution with F = 0,G / 0. 
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Uncoupled Case 

For completeness, it may be worthwhile to write down the solution of the uncoupled field Eq. ©, i.e. 
when d = e = 0. It is easily shown that in this case the solution (|172|) reduces to 

_ / 



a i 
4ci 



1 ± 



m cn 


[v / 2-m( x + 2; o),m_ 




y/(2-m) (yi- m + dn 




) 



(182) 



and one is at the transition temperature T = T c since the parameters oi,6i,ci satisfy b\ = Aa\c\. At 
m = 1, (j) is a constant, which is either or <j) c = yf®*, i.e. at m = 1, the field 4> stays at one of the 
degenerate minima of the potential. 

Energy: Corresponding to the periodic solution (|172|) . while the energy is same, the value of C differs 
depending on if G = 0, F ^ or if F = 0, G ^ 0. For example, if G = 0, F ^ 0, the energy E and the 
constant C are given by 

E = (1 - m){A* + tf 2 )D / 2 " Sn2(y ' m) , 



-21c [VI — m + dn(y, m)] 2 

C = [ai - 6iF + Cl F 2 ]F 2 - \ (A 2 + H 2 )D 2 , (183) 

8 

where y = + ^o)- ON the other hand, if F = 0, G ^ 0, then C is given by 

C = [a 2 - b 2 G + c 2 G 2 ]G 2 - \{A 2 + H 2 )D 2 , (184) 

8 

Unfortunately, we are unable to solve the integral (|183|) analytically and hence are unable to calculate 
the corresponding soliton interaction energy. However, the energy of the corresponding solution [Eq. I|18ip ] 
is of course zero, being a constant solution. 

Solution X 

We now discuss a periodic solution which in the uncoupled limit and even at m = 1 is distinct from 
the well known kink solution (|32j) . This periodic solution is given by 

j ) _ F | Asn[D(x + x ),m] _ Q | Hsn[D(x + x ),m] 

1 + Bsn[D(x + xq), m] ' 1 + Bsn[D(x + xq), m] ' 

provided the following eight coupled equations are satisfied 

2aiF - 3hF 2 + 4ciF 3 + dG 2 + 2eFG 2 = -2BAD 2 , (186) 
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2a x A - SbxFA + 12ciF 2 A + 2dGH + AeFGH + 2eAG 2 = A[6B 2 - (1 + m)]D 2 , (187) 

- 3hA 2 + YlcxFA 2 + dH 2 + 2eFH 2 + AeAGH = 3AB[(1 + m) - 2B 2 ]D 2 , (188) 

2 Cl A 2 + eH 2 = (1 - B 2 )(m - B 2 )D 2 , (189) 

2a 2 G - 3b 2 G 2 + 4c 2 G 3 + 2dFG + 2eGF 2 = -2BHD 2 , (190) 

2a 2 F - 66 2 G# + l2c 2 G 2 H + 2(iF# + 2dAG + AeFGA + 2eHF 2 = H[6B 2 - (1 + m)]D 2 , (191) 

- 3b 2 H 2 + 12c 2 GH 2 + 2cL4# + 2eGA 2 + AeAFH = 3HB[(1 + m) - 2B 2 ]D 2 , (192) 

2c 2 H 2 + eA 2 = {I- B 2 ){m- B 2 )D 2 . (193) 

Six of these equations determine the six unknowns A, B, H, D, F, G while the other two equations, give 
two constraints between the eight parameters ai ;2 , &i,2, c± j2 , d, e. In particular, A and H are given by 

2 _ (l-B 2 )(m-B 2 )D 2 (2c 2 -e) 2 _ (1 - B 2 )(m - B 2 )D 2 (2 Cl - e) 

(4c lC2 -e 2 ) ' (4c lC2 -e 2 ) ' ^> 

Unfortunately, it is not very easy to identify a case when the analysis is somewhat simpler. One possible 
case is when F = 0. In that case, it is easily seen that five of the above equations determine the five 
unknowns A, B, G, H, D while the remaining three equations give three constraints. Further, solutions also 
exist when 2c\ = 2c 2 = e as well as when b 2 = c 2 = 0. However, because of lack of simplicity of these 
solutions, we do not write them here explicitly. 
Uncoupled case at m = 1 

There is one case when the analysis is simple and an explicit solution can be written down. In particular, 
in the uncoupled case, at m = 1 it is easily shown that 



1 + tanh 



t {x + xq 



(1 + B) 44== M, (195) 



4c i \ 1 + B tanh 



' hi 



(x + X ) 



is an exact solution to the above equations. Note that in this case 

b 2 =4a lCl , A = F(l-B), Fb 1 = a 1 (l + B), 2D 2 = oi , (196) 

so that one is again at T = T c . Thus one has generalized the well known kink solution at T = T c by 
obtaining a one parameter family of kink solutions, characterized by the parameter B. In the limit B = 
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this solution goes over to the well known kink solution but is different otherwise. In particular, while as 
x — * ±00, both solutions (|32|) and (|195|) go to 4> c = and respectively, the difference is in how they 
behave at x = 0. In particular, while the kink solution (|32|) goes to 
(|196|) goes to <j> = (1 + ^ n other words, with the new solution, the center of the kink can be 

anywhere between and (j) c , while for the conventional kink solution (|32j). the center of the kink is always 

3/2 

at 4>max- However, the energy of all the kink solutions is identical, i.e. Ekink = and is independent 

of B. 

Energy: Corresponding to the periodic solution (|185|) the energy E and the constant C are given by 
E = (A 2 + H*)D C K ^(y,m)dn 2 (y,m)^ 



tt^, the new kink solution 

2ci ' 



2K [l + Bsn(y,m)Y 



C 



F+ b) 



• ^) • - C 7 ' 77) +e ( G+ f)' 



"2 &2 ( G + |) + C 2 (g + |) 2 + d (F + I) 



(197) 



where y = D(x + xo). Unfortunately, we are unable to solve this integral analytically and hence are 
unable to calculate the corresponding soliton interaction energy. However, the energy of the corresponding 
hyperbolic (bright-bright) soliton solution [Eq. (195)] 



F ^ Atanh[D(x + x )] H taoh[D(x + x )] 



l + Btanh[D(x + x )] 
is easily computed. We find that 

Esoliton ~ 



1 + Btanh[D(x + x ) 
4(A 2 + H 2 )D 



3(1 -£ 2 ) 2 ' 

Note that this solution exists when 2c\ > e, 2c2 > e and Ac\C2 > e 2 . 



(198) 



(199) 



4 Solutions when b\ 7^ 4aiCi (in the Analogous Uncoupled Case): Solu- 
tion XI 



So far, we have considered ten solutions, all of which, in the uncoupled limit correspond to b\ = 4aiCi, i.e. 
T = T c . Now we will discuss a case when the uncoupled limit at m = 1 corresponds to T / T c , i.e. one is 
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either above or below the transition temperature T c and as a consequence one has a pulse lattice solution. 
In particular, it is easily shown that 

Adn[D(x + x ),m] = Hdn[D(x + x ), m] 

9 l + Bdn[D(x + x ),m] ' ^ 1 + Bda[D(x + x ), m] ' 

is an exact solution to the field Eq. (J2J) provided the following eight coupled equations are satisfied 

2a x F - 3hF 2 + 4ciF 3 + dG 2 + 2eFG 2 = 2(1 - m)BAD 2 , (201) 

2aiA - + 12ciF 2 ^ + 2dGH + AeFGH + 2eAG 2 = A[2 - m - 6(1 - m)B 2 ]D 2 , (202) 

- 3hA 2 + 12ciFA 2 + dH 2 + 2eFiJ 2 + AeAGH = -3AB[(2 - m) - 2(1 - m)B 2 ]D 2 , (203) 

2ciA 2 + eH 2 = (B 2 — 1)[1 — (1 — m)B 2 ]D 2 , (204) 

2a 2 G - 3b 2 G 2 + 4c 2 G 3 + 2dFG + 2eGF 2 = 2(1 - m)BHD 2 , (205) 

2a 2 # - 66 2 Gi? + 12c 2 G 2 H + 2dFi? + 2dAG + 4eFGyl + 2ei?F 2 = H[(2 - m) - 6(1 - m)B 2 ]D 2 , (206) 

- 36 2 # 2 + !2c 2 GH 2 + 2(L4F + 2eG^4 2 + AeAFH = -3HB[(2 - m) - 2(1 - m)B 2 ]D 2 , (207) 

2c 2 H 2 + eA 2 = (B 2 — 1)[1 — (1 — m)B 2 ]D 2 . (208) 

Six of these equations determine the six unknowns A, B, H, D, F, G while the other two equations, give 
two constraints between the eight parameters a\ 2, &i 2, ci 2 5 d, e. In particular, A and H are given by 

A 2 = (B 2 -l)[l-(l-m)B 2 }D 2 (2c 2 -e) r2 = (B 2 - 1)[1 - (1 - m)B 2 )D 2 (2 Cl - e) 
(4cic 2 -e 2 ) ' (4cic 2 -e 2 ) 

from where it follows that B 2 must satisfy the bound 

KB 2 < — !— . (210) 
1 — m 

So far, we have not been able to obtain explicit solutions to these equations in an elegant form except at 
m = 1 which we now present. 
m=l 

In case m = 1, the solution (|200|) goes over to the pulse solution 

j ) - Fl Asech[D(x + x ),m] ^ - g \ Hsech l D ( x + x Q ),m] 
1 + Bsech[D(x + xq), m] ' 1 + Bsech[D(x + xo),m] 
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where A, H are now given by 

2 _ (£ 2 -l).P 2 (2c 2 -e) (ff 2 -l)P 2 (2 Cl -e) 
^ " (4c lC2 -e 2 ) ' H ~ (4c lC2 -e 2 ) ' ^ 

Several solutions are possible in this case, which we discuss one by one. 
(i) F = G = 0,m = 1 

In this case one finds that while A,H are given by Eq. (|212[) . the inverse characteristic length D is 
given by 

D 2 = 2oi , 01 = a 2 , (213) 

and B is determined from the relation 

dH 2 + 6AB ai = 3hA 2 . (214) 
Further, the parameters satisfy the constraint 



(2d + 36i)(2c 2 - e) - d(2 Cl - e) = 3W(2ci - e)(2c 2 - e) . (215) 



Interesting case of 6 2 = c 2 = 

In this physically interesting case, one finds that 



A^^-- 1)D \ tf.t^!, (216) 
e e z 

while Eq. ()213|) is still valid. However, B is now given by a simpler expression 

B2 =(2^k)' < 217 > 

and the constraint has a simpler form 

d(e + 2ci) +36ie = 0. (218) 

Special case of 4cic 2 = e 2 

The above solution continues to exist even in the case 4cic 2 = e 2 provided 2c\ = 2c 2 = e. Further, 

while Eq. (|213|) is still valid, A, H are no more given by Eq. (|212|) . but they satisfy the constraint relations 

e(A 2 + H 2 ) = (B 2 — l)D 2 , (2d + 3h)A 2 = dH 2 + 3b 2 HA, 2{dA + 3B ai ) = 3b 2 H , (219) 
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from where one can determine A, H and B. 
(ii) G = 0,F ^0,m = l 

In this case one finds that while A, H are given by Eq. (|212|) . D and F are now given by 



, 36i + J (96? - 32aici) , 
D 2 = 3hF-4 ai , F = <j) c = L . (220) 

OCl 

From the remaining three equations, one can determine B and further one has two constraints between the 
various parameters. This solution is also valid in case 62,02 = or when 2c\ = 2c2 = e with appropriate 
constraints which can be easily worked out. 
(iii) F = 0,G / 0,m = 1 

In this case one finds that the solution exists only if d = 0. Further, A,H are still given by Eq. 1)212(1 . 
while D and G are now given by 



36 2 + J(%1 - 32a 2 c 2 ) 
D 2 = 3b 2 G - 4a 2 , G = i> c = ^— -. (221) 

8c 2 

From the remaining three equations, one can determine B and further one has two constraints between 
the various parameters. This solution is however not valid either in case b 2 , c 2 = or when 2c\ = 2c 2 = e. 
Uncoupled case with m = 1 

Finally, for completeness, we write down the solution in the uncoupled case when m = 1. In this case, 
it is not difficult to show that one of the solution is given by 

J 4sech[v / 2ai(x + xq)] 



1 + -Bsech[y / 2a7(x + xq)] ' 

where 



(222) 



D- . 61 , A= 2Q1 . (223) 
'b\ — Aaicx Jb\-Aa\ci 



The solution (|222j) can also be written in the form |24j 

Tax 



(224) 



b\ + y (b\ — 4oici) cosh[yJ2a\(x + x )\ 
Note that this solution exists provided b\ > Aa\C\, i.e. below the transition temperature T < T c . As 
pointed out at several places in the literature, this solution is a pulse solution around the local minimum 

= 0. 
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The uncoupled equations also have another solution given by 



F Asech[D(x + x ),m] 

l + Bsech[D(x + x ),m] ' 



where 

3bi + J{%\ - 32aici) . 

F = 4> c = L = 36iF - 4ax , 

8ci 

B= Z^" 61 , A = - ; 3M? - 4Q1 ■ (226) 
^(6? _ 2b lCl F) y/(% - 2b lCl F) 

Note that this solution is valid only if b\ < 4oiCi, i.e. above the transition temperature T > T c . This 

follows from the fact that b\ > 2b\C\F so that A, B are real. As pointed out at several places in the 

literature, this solution is a pulse solution around the local minimum <p = <p c = F where 4> c is as given by 

Eq. ©. 

Summarizing, we have obtained eleven solutions to the coupled Eq. © in case a\ > out of which ten 
solutions, for the corresponding uncoupled limit, are at T = T c (i.e. b\ = ka\c\) while one solution is for 
T > Tq or T < T c depending on the value of the parameters. 

5 Solutions when a\ < 

Before concluding this paper, we consider a few solutions of the coupled Eq. (J2J) in case a\ < 0. As 
pointed out in Sec. II, the uncoupled model with a\ < is also a model for first order transition with the 
temperature now being controlled by the parameter b\. We will see that most of the solutions which we 
have obtained in the case a\ > 0, are also valid when a\ < but with suitable modifications of parameters. 
Solution I with F = G = 

In Sec. Ill we saw that with at > 0, solution I as given by Eq. (jlUj) is not possible in case F = G = 0. 
However, it turns out that if a\ < then such a solution is indeed possible. A and B are still given by Eq. 
()19|) while D is now given by 

2 2|oi| 



D 2 = -r^- , (227) 
(1 + m) 



while the three constraints are 



ai = a2 , n h \b x =U\ ^1 = ^4. (228) 

a (2c 2 - e) 

34 



Solution I with G = 0, F / 

Similarly, by looking at Eqs. (|2U|) and (|21|1 it is easy to convince oneself that solution (|1U|) with G = 
and -F = <f) c \ ox F = <j) C 2 exists in case a\ < and the parameters satisfy essentially similar relations to 
Eqs. (|19|) to (|21|) with appropriate modifications. Further, such solutions continue to exist in the special 
case of e 2 = 4ciC2- However, such solutions do not exist in the physically interesting limit of 62 = C2 = 0. 
Note that <p C 2 are as defined by Eqs. © and ©. 

Solution I with F = 0, G / 

Similarly, by looking at Eqs. (|28|) and 1|29|) it is easy to convince oneself that solution with F = 
and G = ^ m exists in case 01 < and the parameters satisfy essentially similar relations to Eqs. (|27j) to 
(|29|) with appropriate modifications. But, such a solution does not exist either when e 2 = 4ciC2 or when 
62 = C2 = 0. Note that ip m ,ip c are essentially same as (j) m ,(pc but with the replacement of ai,b\,ci by 

Solution II with F = G = 

In Sec. Ill we had seen that solution 1|38|) with F = G = and m < 1/2 is not possible in case a\ > 0. 
It is however easy to see that if a\ < 0, then such a solution is clearly possible (with m < 1/2) provided 
relations (|47|) to ()49|) are satisfied. However, like other solutions of type II (and also III) with a\ > 0, no 
solutions of type II with a\ < are possible in case either e 2 = 4ciC2 or if 62 = C2 = 0. 

Solution II with G = 0, F / 

Similarly, by looking at Eqs. ()5U|) and (|51|) it is easy to convince oneself that solution l)38|) with m > 1/2, 
G = and F = <p c i or F = (j) C 2 exists in case a% < and the parameters satisfy essentially similar relations 
to Eqs. I|47|). (fST)|) and ()51j) with appropriate modifications. 

Solution II with F = 0, G / 

Similarly, by looking at Eqs. ()56|) and (|57|) it is easy to convince oneself that solution l)38|) with m < 1/2, 
F = d = and G = exists in case a\ < and the parameters satisfy essentially similar relations to 
Eqs. (|56|). (|57|) and (|59|) with appropriate modifications. 

Solution III with G = 0, F / 

Looking at Eqs. ((77)) and ifTS")) it is easy to convince oneself that solution ^tyi with G = and F = 4> c \ 
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or F = 4> C 2 exists for any m (0 < m < 1), in case a\ < and the parameters satisfy essentially similar 
relations to Eqs. (|74jl. (|77jl and 1|78|) with appropriate modifications. 
Solution VIII with F = G = 

In Sec. Ill we observed that with a\ > 0, solution VIII as given by Eq. (|151j) is not possible in case 
F = G = 0. However, it turns out that if a% < then such a solution is indeed possible. A and H are still 
given by Eq. (|161j) while the inverse characteristic length D is now given by 

D 2 = t (229 ) 

(2-m) v ; 

while the three constraints are 



ai = a2 , 21b\b x =U\ ^1 = 1^-4. (230) 

a (2c 2 - e) 

Solution VIII with F = 0, G ^ 

While no solutions are possible for G = 0, -F 7^ in case ai < 0, it turns out that solutions are indeed 
possible in case F = 0, G 7^ 0, a\ < provided d = and 5 = 1. In this case A, H are given by Eq. (|161|) 
while G and D are given by Eq. (|165|) . Further there are three constraints given by 

ZbiA = AeGH , 3H 2 (b 2 - 4c 2 G) = 2eGA 2 , 4(|ai - eG 2 ) = (2 - m)D 2 . (231) 

Solution IX with F = G = 

In Sec. Ill we also observed that with a± > 0, solution IX as given by Eq. (|172|) is not possible in case 
F = G = 0. However, it turns out that if a\ < then such a solution is indeed possible. A and H are still 
given by Eq. (|161j) while D and the three constraints are again given by Eqs. (|229|) and (|23fl|) . respectively. 

Solution IX with F = 0, G / 

Since solutions VIII and IX satisfy same field equations, it is then clear that a solution with F = 0,G / 
0, a\ < will also exist in this case. 
Solution XI 

In the last section we obtained a pulse-like periodic solution in case a\ > 0. It is amusing to note that 
the same solution continues to be valid even when a\ < 0. In particular, note that in the uncoupled limit 
and at m = 1 the solutions for a\ < are again pulse solutions around the local minimum. Specifically, if 
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b\ > (<) 0, then the pulse solution is around cj) C 2 (<fi c i)- in both cases, the solution is given by Eq. ()225j) 
but with F now being either <\> c \ or <\> c i where (/)&, (p C 2 are as given by Eqs. (J5J) and ©, respectively, and the 
values of other parameters like D, B, A, H are again given by Eq. (|226|) but with appropriate modification. 

6 Conclusion 

We have systematically provided an exhaustive set of exact periodic domain wall solutions for a coupled 
asymmetric double well model (with and) without an external field. We found ten solutions at the transition 
temperature (T = T c ) and one solution each for T > T c and T < T c . For the (physically interesting) special 
case when hi = C2 = there is no nonlinearity in the tp field. However, due to the coupling with the <f> field 
which has explicit nonlinearity, an effective nonlinearity is induced in the tp field thus enabling soliton-like 
solutions. When we set the coupling parameters d = e = the soliton solutions cease to exist in the ip 
field, as expected (for b<i = C2 = 0). We emphasize that the solutions found here are quite different from 
the ones we found for the first order transition in the coupled <ft e model |22| or the second order transition 
in coupled symmetric double well (</> 4 ) model [5T] . 

It would be instructive to explore whether the different solutions reported here in Sec. Ill are completely 
disjoint or if there are any possible bifurcations linking them via, for instance, analytical continuation. We 
have not tried to carry out an explicit stability analysis of various periodic solutions. However, the energy 
calculations and interaction energy between solitons (for m ~ 1) could provide useful insight in this 
direction. Unlike the coupled 4 and cases |21l 122 j . we have not yet succeeded in finding a solvable 
discrete analog of an asymmetric double well model. 

Our results are relevant for a wide range of physical situations including structural transformations 
(with coupling of strain and shuffle modes) [HI El j liquid crystals [H] , hydrogen bonded chains |1J 03 H3 [7] 
and field theoretic contexts [S]. For structural phase transformations the solutions provide novel domain 
wall arrays, e.g. periodic antiphase boundaries and twin boundaries. Similarly, for hydrogen bonded chains 
these solutions represent new periodic nonlinear excitations. 
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